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ABSTRACT 


Satellite  data  first  became  available  in  the  1960's 
making  solar  budget  calculations  more  feasible.   In  this 
study,  satellite  measurements  of  reflected  solar  insolation 
from  NIMBUS  II  were  combined  with  simultaneous  pyrheliometer 
readings  of  transmissivlty  over  North  America  to  provide  the 
basis  for  computation  of  the  atmospheric  absorptivity.   In 
the  statistical  parameterization  of  the  absorptivity,  such 
readily  available  data  as  water  vapor  mass  and  satellite  re- 
flectance were  introduced  as  independent  variables  in  the 
regression  equations.   Correlation  coefficients  between  obser- 
ved regression-determined  absorptivities  and  transmisslvities 
were  computed  based  upon  an  initial  total  of  235  data-samples, 
and  the  results  interpreted.   The  conclusion  was  twofold: 
(1)  A  new  suggested  insolation-budget  model  is  proposed  for 
computation  of  absorptivity;  (2)  the  desirability  of  deriv- 
ing a  statistical  parameterization  for  cloud-cover  at  grid 
points  in  terms  of  satellite-observed  reflectances  and  window- 
channel  temperatures  is  shown. 
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I.   INTRODUCTION 

The  solar  energy  cycle  Is  of  vast  importance  to  the 
meteorological  community  since  addition  of  heat  is  a  pri- 
mary source  of  energy  in  atmospheric  circulations.   There- 
fore, a  proper  accounting  of  absorptivity  and  transmissiv- 
ity within  the  atmosphere  is  vital  to  improved  modeling 
and  numerical  forecast  procedures. 

This  study  is  designed  to  correlate  satellite  measure- 
ments of  reflected  solar  insolation  from  the  Medium  Resolu- 
tion Infrared  Radiometer  (I-iRIR)  on  NIMBUS  II  and  simultan- 
eous pyrhellometer  readings  of  transmissivity  over  North 
America.   This  research  was  directed  towards  the  parameter- 
ization of  the  atmospheric  absorptivity  and  transmissivity 
in  terms  of  available  satellite  albedo  observations.   In 
addition  such  other  readily  available  data  such  as  water 
vapor  mass  and  cloud-cover  may  be  introduced  for  perfect- 
ing the  parameterization. 

The  limits  of  the  study  were  confined  to  the  United 
States  and  Canada  where  data  were  readily  available. 

Studies  of  the  disposition  of  solar  Insolation  prior  to 
the  acquisition  of  satellite  data  have  been  reported  by 
Gabites  (1950),  London  (1957)  and  Sellers  (1965).   Each  of 
these  researchers  gave  estimates  of  global  albedo  of  close 
to  35;t.   Vonder  Haar  (1968)  using  satellite  data,  has  shown 
the  global  albedo  to  be  close  to  2>Q%.      In  view  of  this 
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reduction  of  the  albedo  of  %   from  the  previous  estimates, 
one  should  look  for  revised  transmisslvity  and  absorptivity 
parameterizations  which  would  account  for  the  insolation 
balance . 

It  would  seem  that  successful  parameterization  of  ab- 
sorptivity and  transmisslvity  could  be  accomplished  over 
data-rich  areas  such  as  existed  over  North  America  during 
the  NIMBUS  II  period  of  this  study.   In  the  long-term,  one 
would  hope  to  describe  these  variables  of  the  solar  budget 
in  terms  of  readily  accessible  satellite  observations,  and 
ultimately  expand  the  parameterization  to  oceanic  and  other 
data-sparse  areas. 
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II .   DATA  SOURCES  AND  PROBLEM  OUTLINE 

A.  GENERAL 

The  data  used  in  this  study  was  from  three  basic  sources 
and  covered  a  five-day  period  from  15  July "1966  to  20  July 
1966.   The  surface  data  were  compiled  at  the  ^7  pyrhelio- 
metric  stations  shown  in  Fig.  1.   The  locations  were  chosen 
to  provide  a  broad  coverage  of  solar  transmission  effects 
on  the  North  American  continent  within  the  limits  25°  north 
latitude  to  59°  north  latitude.   Several  other  types  of 
data  compilations  (described  below)  were  required  to  perform 
computations  of  solar  reflectance,  transmittance,  and  ab- 
sorptivity at  pyrheliometric  stations  at  the  effective  time 
of  the  NIMBUS  II  scan  time.   Thus  all  observational  data 
had  to  be  interpolated  to  local  solar  times  corresponding 
to  the  effective  scan  time  of  the  NIMBUS  II  satellite  at 
the  pyrheliometric  station.   Solar  time  is  computed  every- 
where on  the  basis  that  1200  solar  time  corresponds  to  the 
highest  solar  elevation  of  the  day. 

B.  NIMBUS  II  REFLECTANCE  DATA 

NIMBUS  II  solar  reflectance  values  were  available  for 
this  study  only  in  hard  copy  printout  format.  The  NIMBUS 
II  satellite  had  a  nearly  sun-synchronous  polar  orbit  with 
the  ascending  node  at  11:28  (equatorial  solar  time).   The 
satellite  scanned  a  swath  perpendicular  (Fig.  2)  to  its 
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Figure  1  -  Mercator  map-projection  of  data  region  showing 
pyrheliometer  station  locations  used  in  this  study.   De- 
tailed latitude  and  longitude  values  for  each  station  are 
located  in  Appendix  Table  A. 
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orbit,  and  extending  between  the  horizons.   Each  scan  spot 
is  approximately  a  30-mlle  square  at  the  subsatellite  point 
and  increases  in  size  with  increasing  nadir  angle  of  the 
scan  spot  relative  to  the  satellite. 

The  NIMBUS  II  Medium  Resolution  Infrared  Radiometer 
(MRIR)  was  designed  to  scan  in  the  five  wavelength  regions 
(listed  below)  with  each  region  accompanied  here  by  a 
brief  description  of  the  observational  purpose  of  channel 
observation  (e.g.,  NIMBUS  II  User's  Guide,  1966). 

6.4  to  6.9  microns  -  This  channel  covers  the  6.7 
micron  water  vapor  absorption  band.   Its  purpose  is 
to  provide  information  on  water  vapor  distribution 
in  the  upper  troposphere  and,  in  conjunction  with 
the  other  channels  to  provide  data  concerning  rela- 
tive humidities  at  these  altitudes. 

10  to  11  microns  -  Operating  in  an  atmospheric 
"window,"  this  channel  measures  surface  or  near- 
surface  temperatures  over  clear  portions  of  the  at- 
mosphere.  It  also  provides  cloud  cover  .and  cloud 
height  information  (day  and  night). 

14  to  16  microns  -  This  channel,  centered  about 
the  strong  absorption  band  of  CO2  at  15  microns, 
measures  radiation  which  emanates  primarily  from  the 
stratosphere.   The  information  gained  here  is  of  pri- 
mary importance  in  following  the  seasonal  stratos- 
pheric temperature  changes. 
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Figure  2  -  Mercator  map  showing  NIMBUS  II  orbital  paths 
as  indicated  by  the  solid  lines  — .   The  limits  of  the 
swath  coverage  are  indicated  by  the  dashed  lines  — . 
Overlapped  regions  are  indicated  by  the  cross  hatching//. 
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5  to  30  microns  -  This  channel  measures  the  emit- 
ted long  wavelength  infrared  flux  and,  in  conjunction 
with  the  reflected  solar  radiation  channel  furnishes 
data  on  the  heat  budget  of  the  planet. 

0.2  to  4.0  microns  -  This  channel  covers  more  than 
99$  of  the  solar  spectrum  and  yields  information  on 
the  intensity  of  reflected  solar  energy  from  the  earth 
and  its  atmosphere. 

These  five  channels  are  frequently  termed  channels  1, 
2,  ....  5,  respectively,  of  the  MRIR  experiment  of  NIMBUS 
II.   Of  the  various  channel-data  available,  only  channel  5 
relating  to  solar  reflectance  from  the  earth  and  its  at- 
mosphere was  utilized. 

It  would  perhaps  have  been  a  useful  indirect  measure 
of  cloud  cover  to  have  used  the  "window"  channel-equivalent 
blackbody  temperature  data  for  the  same  observation  times. 
However,  synoptic  surface  temperatures  at  the  pyrhellometric 
stations  would  have  been  necessary  to  give  a  concise  in- 
dication of  the  cloud  amount  and/or  height,  so  that  the 
channel  2  data  were  not  analyzed  here. 

The  effective  radiant  reflectance  R  in  the  direction 
of  the  satellite  from  the  reflecting  surface  filling  the 
field  of  view  of  the  radiometer  and  illuminated  by  unatten- 
uated  solar  radiation  is  defined  (p.  57,  NIMBUS  II  Users' 
Guide,  1966)  by 
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N  I1 .J' ,Z' 
BII'.J'.Z')  =  ^ ' *-  (1) 

Sr*  cos  Z' 

m 

r2 


where  I*,J'  are  the  coordinates  of  the  scan  spot,  Z'  is  the 
zenith  angle  of  the  sun  and 

H  =  ^i  (2, 

r^ 

is  the  date-adjusted  value  of  the  solar  constant  S  (correc- 
ted to  17  July  1966).   The  parameter  N  in  (1)  is  the  solar 

— ?     — l 
intensity  or  radiance  in  watts  m   ster   reflected  back 

to  the  satellite,  and  is  here  assumed  independent  of  the 
viewing  angles  (nadir  and  azimuth)  of  the  satellite  rela- 
tive to  the  scan  spot.   The  latter  assumption  is  not  ac- 
tually valid  as  has  been  shown  by  Raschke  et  al.  (1973) 
who  showed  the  distribution  of  the  bidirectional  reflec- 
tance over  typical  land  spots  to  be  dependent  upon  both 
angles.   However,  our  scan-spot  reflectances  were  computed 
by  NASA  according  to  Eq.  1.   The  hard-copy  reflectance  data 
obtained  from  NASA  were  in  the  form  of  averages  composited 
to  gridpoints  on  a  Mercator  map  projection  of  scale  1:10M 
(Fig.  3)  and  with  a  grid-mesh  (I, J).   This  (I, J)  grid-mesh 
may  be  formulated  in  terms  of  longitude  X  and  latitude  <P 
by  means  of  the  transformation: 

I  -  1  =  0.8(150°  -  X°.)  (3) 
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Figure  3  -  Mercator  map-projection  of  data  region  used  in 
analysis  of  data  considered  in  this  study.   Equally  spaced 
I-  and  J-  intervals  are  indicated  along  the,  bottom  and 
left  boundaries  of  the  rectangular  area.   Corresponding 
longitude  and  latitudes  are  shown  along  the  opposite  boun- 
daries. 
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j  -  1  =  o.QilM)   m  (1  +  slnj>)  (4) 

71  COS  0    '  N  ' 


Some  of  the  NIMBUS  II  data  orbits  used  on  16  July  1966  are 
sketched  in  Fig.  2. 

Within  any  grid  space  centered  at  (I, J)  the  number  of 
scan  spots  may  vary  from  five  to  ten  depending  on  the  nadir 
angle  to  the  scan  spot,  the  actual  earth-size  of  the  grid- 
area  (which  is  determined  by  the  map  factor  m  =  sec  4>  ,  dg  = 
dmap/m) ,  and  the  extent  of  swath  overlap  at  the  scan  spot 
between  successive  orbits. 

The  mean  reflectance  value  obtained  at  each  grid  point 
(I, J)  is  therefore  an  average  composited  from  all  of  the 
individual  scan-spot  reflectances  R(I, J, I ' , J1 ,Z' ) .   Here 
I*,J'  are  individual  scan  spots  contained  in  the  larger 
grid-mesh  I, J  of  the  map  (Fig.  3)  within  the  grid  square. 
The  averaged  value  R(I,J,Z)  is  the  mean  solar  reflectance 
to  space,  computed  after  averaging  over  all  the  scan  spots 
I', J1  within  I, J  and  is  a  function  only  of  the  gridpoint 
identifiers  (I, J)  and  of  Z,  the  mean  solar  zenith  angle  at 
(I, J).   R(I,J,Z)  as  thus  determined  is  termed  the  direc- 
tional reflectance  at  the  grid  point,  the  direction  being 
that  of  the  mean  solar  zenith  angle  to  the  center  of  the 
grid  square.   Thus  after  the  averaging  reflectances  to  grid 
points  has  been  accomplished,  the  resulting  reflectance 
R(I,J,Z)  is  given  by: 

R(I.J.Z)  =  rriT/  I  S/r  \2  cos  zl  (5) 


""/    [Sfol)2  C°S  Zl 


19 


where 

5  =  solar  constant  =  2.0  langleys  mln"1 

cos  Z  =  mean  cosine  of  the  solar  zenity  angle  sub- 
tended by  point  I, J  at  station. 
cos  Z  =  cos  <t>   cos  6  cos  h  +  sin  <f>   sin  6  (6) 

<t>  =  station  latitude 

6  =  solar  declination 

h     =  hour  angle  at  effective  time  of  satellite 

scan. 
N     =  mean  reflected  solar  radiance  from  (I, J)  to 

the  satellite. 
(rm/r)=  ratio  of  the  mean  distance  to  the  sun  and  the 

distance  to  the  sun  at  the  required  time. 

1.   Effective  Tine  Calculation 

The  effective  Greenwich  wean  Time  of  the  satellite 
reflectance  observation  at  any  station  was  computed  using 
a  linear  interpolation  scheme  between  the  times  along  the 
swath  over  the  station  and  orthogonal  to  successive  satel- 
lite tracks,  and  the  stations  distances  relative  to  the 
tracks  as  shown  by  Fig.  *K   Thus  Eq.  (7)  is  used  to  inter- 
polate for  effective  time  over  the  station  as  shown  in 
Fig.  *K 

Effective  scan-time  =  (^)(T2)  +  (2i=4£)  (t  )      (7) 

If  the  station  was  in  a  region  scanned  only  by  one  swath, 
the  time  I,  was  read  directly  off  the  nearest  track  in 
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Figure  k.      Interpolated  Scan-Time 
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Fig.  2.   The  effective  time  T  was  then  converted  to  true 
solar  time  for  the  station  by  the  following  equation: 

S-T-  - I(GMI)  -  iJtSS  *  (aRS>  <8) 

Note  the  hour  angle  h  depends  upon  the  solar  time  of  the 
scan  as  given  by  Eq.  9  and  has  the  form: 

h  =  S.T. (hours)  -  12  *  15°  (9) 

2.   Station  Reflectance  Values 

The  pyrhellometer  station  location  (Fig.  1)  falls 
within  groups  of  four  grid  points  shown  in  Fig.  5*   There- 
fore, a  linear  Interpolation  scheme  was  used  to  determine  a 
value  for  the  exact  station  location  (I+Al ,    J+AJ)   located 
between  the  four  adjacent  grid  points  at  which  the  reflect- 
ance values  H(A),  R(B),  R(C),  R(D)  were  available.   The 
reflectance  value  at  the  location  I+Al ,  J+4J  (the  station 
location)  is  given  by  means  of  the  standard  spatial  inter- 
polation formula  . 

R  =  (1-AJ)    [(1-4I)R(A)+4I*R(3)]  +4J  [ (1-4I)R(D) 

+4I*R(C)]  (10) 

Here  R(A)  is  the  reflectance  at  point  A,  R(B)  is  that  at 

B,  etc. ,  as  indicated  by  the  grid  array  in  Fig.  5« 

C.  SURFACE  TRANSMITTANCE  OF  SOLAR  INSOLATION 

Surface  transmit tance  of  solar  radiation  was  based  on 
pyrhellometer  readings  recorded  at  the  stations  in  Fig.  1. 
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Figure  5-   Reflectance-Interpolation  to  Stations, 
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These  records  provided  on  magnetic  tape  are  listings  of  the 
hourly  accumulations  of  solar  Insolation  observed  at  the 
end  of  each  true  solar  hour.   The  taped  record  was  kindly- 
furnished  with  the  technical  assistance  of  the  Commander, 
Naval  Weather  Service.   The  values  listed  were  pyrhelio- 
meter  readings  of  the  total  hemispheric  insolation,  with 
the  pyrheliometer  calibrated  to  conform  to  a  solar  constant 
of  2.0  langleys  min~l. 

The  locations  of  pyrheliometrlc  stations  were  converted 
from  latitude  0  and  longitude  X  to  the  (I, J)  coordinate 
system  in  accordance  with  Eqs.  3  and  4.   The  mean  station 
transmlttances  were  determined  at  the  effective  time  of  the 
satellite  scan.   This  was  done  to  give  a  transmlttance 
coincident  with  the  "simultaneous"  satellite  mean  reflect- 
ance observations  over  this  station.   The  actual  station 
insolation  was  considered  to  be  time-centered  at  0930  solar 
hours  when  it  was  reported  for  the  solar  hour  ending  at 
1000,  etc.   Thus,  it  was  necessary  to  interpolate  an  in- 
solation value  at  the  station  coincident  with  the  effective 
time  of  the  satellite-scan.   This  was  accomplished  by  using 
a  cubic  Lagrangian  time-interpolation  scheme  described 
below. 

The  four  successive  values  of  the  hourly  insolation  that 
bracketed  the  effective  scan  time  were  used  in  the  following 
equation: 


2k 


Q-  =  (T-T2)«(T-T3)«(T-rM^l    (T-T1)»(T-T3)»(T-T*0«Q2 
^    (Tl-T2)*(T1-T3)*(ri-1'4)  +  (l-2-Tl)«(T2-l'3)*(T2-l^J" 

+  (T-Tl)MT-T2)*(T-lM*?3  .  (T-Tl )» (T-I2 )* (T-T3 )«Q>  rvn 
(T3-T1)*(T3-T2)*(T3-t4J  +  1t4-T1)*(T4-'I-2)*(T^-T3)  U1) 

Here 

QG  =  station  insolation  interpolated  to  the 
scan  time. 
T  =  effective  time  of  satellite-scan. 
T1,T2,T3»T^  =  successive  times  of  pyrheliometric 

readings  at  half-hour  solar  clock  times. 
Q1,Q2,Q3,Q^  =  station  insolations  at  the  four  suc- 
cessive times  T1,T2,T3,T4. 

The  Langranglan  fit  was  used  to  reduce  the  effect  of 
transistory  clouds  over  a  short  time  frame. 

D.   PRECIPITABLE  WATER  VAPOR 

Precipitable  water  vapor  was  computed  from  radiosonde 
soundings  nearest  to  the  pyrheliometric  stations.   The 
soundings  were  taken  at  0000GMT  and  1200GMT  starting  at 
0000GI4T  16  July  1966,  and  continuing  through  0000GHT  20 
July  1966.   16  levels  were  employed  using  the  surface  as 
the  bottom  level  and  the  15  standard  levels  at  whole  mul- 
tiples of  50  nib  above  the  surface.   Temperatures  and  rela- 
tive humidities  (RH)  at  these  levels  were  extracted  from 
the  soundings  and  used  to  compute  water  vapor  mixing  ratios 
at  the  indicated  levels  in  the  vertical.   The  computation 
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of  q  was  facilitated,  with  RH  known,  by  use  of  the  Goff- 
Gratch  (cf.,  List,  1958.  p.  350)  formulation  of  the  satura- 
tion vapor  pressure 

log1()ew  =  -7. 90298 (Ts/T-l)+5. 02809  log1Q(Ts/T)- 
1.3816  A   lO-^doH^^d-T/Ts)^)  + 

8.1328  *  10-3  *  (io-3-/+91^9(Ts/T-l).1)  + 

1Og10ews  (12) 

where 

ew  =  saturation  vapor  pressure  over  a  plane  surface  of 
pure  ordinary  liquid  water. 
T  =  absolute  temperature  (°K) . 
Ts  =  steam  point  temperature  (373.l6°K). 
eWs  =  saturation  pressure  of  pure  liquid  water  at  steam 
point  temperature  (1  standard  atmosphere  = 
1013.246  mb). 

Following  the  use  of  Eq.  (12),  the  value  of  q  is  given  by: 

q  =  621.97  ^  =  62l.97(EH)ew  (13) 


p-e 

p-(RH)e 


w 


Then,  precipltable  water  vapor  is  found  from  the  mixing 

ratio  integration 
r300 

„ .  imf    qdp 

&  SIC 
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For  the  times  on  16  July  1966  through  19  July  1966  a  linear 
interpolation  scheme  between  the  1200  and  0000GHT  soundings 
was  used  to  establish  the  precipltable  water  vapor  at  the 
effective  time  of  satellite  scan  at  each  station. 

More  specific  examples  of  the  calculations  performed 
in  order  to  relate  the  satellite  reflectance,  and  simul- 
taneous station-transmittance  through  a  physical  relation- 
ship are  detailed  in  Section  III. 
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III.   DATA  PROCESSING 

A.   THE  SOLAR-INSOLATION  MODEL 

In  Section  II,  a  general  description  was  given  of  the 
data-collection  problem  at  each  pyrheliometric  station  at 
the  effective  time  of  the  satellite  scan.   The  immediate 
objective  of  the  data-analysis  corresponding  to  each  station- 
observation  itfas  to  solve  the  solar-insolation  budget  equa- 
tion for  the  atmospheric  absorptivity. 

FIGURE  6 


Q 


4 


QR|   Top  of 


atmosphere 


Qa 


%| 


aG*%t 


EARTH  ST3HF5CE 

Figure  6.   The  solar-insolation  model 
over  a  pyrheliometric  station.   See 
Eq.  15  and  following  for  explanation 
of  symbols. 


Fig.  6  shows  an  insolation  model  after  Fritz  et  al.  (196^) 
and  Hanson  et  al.  (1967)  which  illustrates  the  terms  involved 
in  the  solar-insolation  budget  at  a  satellite  scan-time  for 
a  given  pyrheliometric  station.   These  terms  are: 


Qs  =  S(rm/r)   cos  Z 


(15) 
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Qs  =  incident  extraterrestrial  insolation  based  upon  a 

solar  constant  S  of  2.00  langleys  min"1  (Qs  is  the 
denominator  of  Eq .  5)« 

Qq  =  pyrheliometrlc  value  of  the  insolation  at  the  sta- 
tion interpolated  to  the  true  solar  time  (Eq.  8) 
of  the  satellite  scan,  according  to  Eq.  11. 

Qr  =  RQs,  which  is  the  reflected  insolation  at  the  sta- 
tion computed  using  the  space- interpolated  reflect- 
ance R(I+4I,  J+4J)  of  Eq.  10. 

Q^  =  absorbed  insolation  in  the  atmospheric  column  above 
the  station. 

aG  =  surface  albedo  of  the  earth  In  the  vicinity  of  the 
station. 

The  midsummer  values  of  ac   were  taken  from  a  North 
American  atlas  reproduced  here  as  Fig.  7  (prepared  by  Kung 
et  al.,  196^).   The  term  1-«g  is  the  fractional  absorption 
of  the  incident  insolation  Qq,  at  earth.   At  each  of  the  k? 
pyrheliometer  stations  used,  the  atlas  value  of  aG  has  been 
listed  next  to  the  station  coordinates  (Appendix  A). 

The  insolation-budget  equation  then  according  to  Hanson 
et  al.  (1967)  becomes: 

Qs  =  QR  +  QA  +  Qo  U-«g)  (l6) 

If  each  term  of  (16)  is  divided  by  Qs,  the  equation  may  be 
rewritten  as 
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Figure  7  -  Summer  surface  albedo  map  for  North  America 
(after  Kung  et  al. ,  1964) 
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qA  =  1  -  qR  -  qG  (1-  <*Q)  (17) 

where : 

qA  =  atmospheric  absorptivity. 

qR  =  Qg/Q-,  is  the  computed  reflectance  using  the  NIMBUS 

II  data-maps. 
qG  =  atmospheric  transmissivity  Qq/Q.5  as  inferred  from 

the  pyrheliometer  reading. 

The  parameter  a„   was  taken  from  the  summer  mean  reflect- 
ance maps  of  Kung  et  al.  (1964)  for  each  station,  and  assum- 
ed constant  during  the  NIMBUS  II  period  15-19  July  1966,  al- 
though a   varied  from  station  to  station  in  accordance  with 
G 

varying  terrain  and  vegetation  characteristics.   The  use  of 
the  aG  in  connection  with  Eq.  1?  was  considered  valid  by 
Hanson  et  al.  (1967)  both  in  clear  and  partly  cloudy  sky 
conditions,  and  this  same  approximation  is  made  here  in 
reference  to  the  computation  of  q^  in  conditions  of  mixed 
cloudiness.   An  attempt  is  made  in  the  subsequent  regression 
analysis  to  isolate  that  part  of  q^  which  depends  upon 
cloudiness.   Statistical  specification  of  q^  similar  to  that 
of  Hanson  will  be  performed  here  using  the  simultaneously 
observed  and/or  calculated  values  of  q£,  qG ,  and  a      during 
the  period  15-19  July'1966,  with  qR  measured  by  NliYiBUS  II 
satellite. 

B.   DATA  TABULATIONS 

For  present  tabulation  purposes,  the  47  insolation  sta- 
tions were  assigned  arbitrary  station  numbers  1,  2,  ...,  47 
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TABLE  I.   TYPICAL  SET  OF  INSOLATIONAL 
AND  RELATED  DATA  AT  PYRHELIOMETER  STATIONS  AT 
THE  EFFECTIVE  TIME  OF  THE  NIMBUS  II  SCAN  ON  17  JULY  1966 


Extra- 

terrest 
Station   Qg, 

No.   lynln"! 


F 
F 


1 
2 

3 

4 

5 
6 

7 

8 

9 
10 

11 

12 

13 

14 

15 

16 

17 

18 

19 
20 
21 
22 
23 
24 

25 
26 

27 

28 

29 
30 

31 
32 

33 
34 

35 
36 
37 

38 

39 
40 
41 
42 

43 
44 

45 
46 

47 


1.7790 
1.8857 

1.9189 
1.7939 

1.8884 

1.8381 
1.7782 

1.7949 
1.6741 
1.7562 

1.8971 
1.8742 
1.8047 
1.8233 
1.8660 
1.8013 
I.6696 
1.7226 
1.7201 
1.6415 

1.7945 
1.8787 
1.7861 
1.7012 
1.7691 
1.7277 
1.5645 

1.4358 
1.6297 
1.5447 

1.6977 
1.5499 
1.7278 

1.7287 
1.6254 

1.7199 
1.6432 
1.6118 
1.6971 
1.6605 

1.6035 
1.6532 
1.629? 
1.7255 
1.7255 
I.636I 
1.6160 


Re- 
flected 

^R 

0.6398 
0.2711 
0.3262 
0.3243 
0.7836 
0.6659 
0.3921 
0.3647 
0.2470 

0.3273 
0.4945 
0.4056 
0.3676 
0.228? 
0.3730 
0.2279 
0.5659 
0.3008 
0.2449 

0.5903 
0.3556 
1.0111 
0.3806 
0.2461 
0.3044 
0.3217 
0.5559 
0.7002 
0.5130 
0.4654 
0.2495 
0.4177 
0.2937 
0.2173 
O.5605 

0.2183 
0.2796 
0.6234 
0.3342 
0.1193 
0.3091 
0.5877 
0.2803 
0.2697 
0.2697 
0.4925 
0.6268 


Pyrhel- 
lometer 


Absorbed  PJJter   Zenith 
(Air)    Vapor    Path 
Qa    U.gin  cm"2   SEC  Z 


0.8319 
1.3564 
1.6974 
1.4725 
0.6680 
0.6116 
1.0544 
1.1184 
1.2202 
1.1706 
1.1978 
1.6545 
1.3676 

1.3959 
1.2186 

1.3395 
-0.0474 
1.1724 
1.3934 
1.3108 
1.4354 
0.2224 
1.2782 

1.0973 
1.2806 

1.3567 
0.3864 

1.0219 
1.1258 
1.0314 
1.2607 
1.0461 
1.2433 
1.2437 
0.4457 

1.1793 
1.2288 

0.7737 
1.1203 
I.6500 
I.6500 

1.7551 
1.1822 

1.2719 

1.2948 
0.9245 
0.6268 


0.4404 

0.4210 
0.1500 
0.2328 
O.5236 
0.6461 
0.5004 
0.5130 
0.3778 
0.4573 
0.4324 
0.1782 
0.3294 
0.4080 
0.5424 
0.4348 
1.1439 
0.4487 
0.3048 

-0.0499 
0.1901 
0.6764 

0.3191 
0.5225 
0.3762 
0.2528 
0.6763 
-0.1330 
0.1598 
0.202? 
0.3766 

0.2535 
0.3648 
0.4666 
0.6861 
0.4992 

0.3315 
0.3308 
0.4218 

0.1553 
-0.1081 
0.4311 
0.3445 
0.3620 
0.3423 
0.3671 
0.4459 


4.8232 

4.5175 
4.3973 
4.4853 
4.3882 
4.3088 
3. 0089 

2.6633 

1.8045 

2.2867 
2.4511 
2.1694 
1.3708 
2.18?8 
3.2688 
1.3032 
3.0095 
1.5689 
1.6416 

1.6987 
I.6033 

4.5055 
2.4676 
1.4247 
1.5437 
3.1032 
2.2638 
2.1568 
1.6994 

2.2699 
1.8045 

2.1985 
1.4729 
1.3883 
3.0243 
1.9092 
2.1046 

2.8319 
2.1268 
1.4620 
1.5630 

1.5231 
1.6260 
1.4731 
1.4731 
1.6994 
3.1851 


1.0883 

1.0267 
1.0089 
1.0792 
1.0252 

1.0533 
1.0887 
I.0786 

1.1565 
1.1024 
1.0205 
1.0330 

1.0727 
1.0618 
1.0376 
1.0748 
1.1596 
1.1238 

1.1255 
1.1794 
1.0?89 
1.0305 

1.0839 
1.1380 
1.0943 
1.1206 

1.2375 
1.3484 

1.1879 
1.2533 
1.1404 
1.2491 
1.1204 

1.1199 
1.1910 
1.1256 
1.1781 
1.2012 
1.1408 

1.1659 
1.2074 

1.1711 

1.1879 
1.1220 
1.1220 

1.1833 
1.1966 
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(see  Appendix  A,  for  station  names,  code  and  geographical 
coordinates).   Thus,  for  each  station  of  Fig.  1,  the  values 
Qs ,  Qq ,  QR  were  computed  using  Eqs.  15.  11.  10;  and  QA 
was  computed  as  the  residual  from  Eq.  16.   A  typical  set 
of  insolation  values  obtained  for  each  of  the  ^7  stations 
on  17  July  1966  is  listed  in  Table  I. 

The  primary  objective  here  at  the  outset  was  to  obtain 
a  statistical  regression  for  the  absorptivity  q^  against 
an  effective  measure  of  the  water  vapor  absorbing  mass  in 
the  air  column  above  the  pyrheliometer  station.   For  this 
purpose  simultaneous  values  of  precipitable  water  vapor  U 
and  of  sec  Z  were  also  listed  in  Table  I. 

The  complete  set  of  raw  data  to  be  used  for  this  re- 
gression analysis  is  tabulated  in  Appendix  B.   There  were 
in  all  ^7  X  5  =  235  simultaneous  values  of  the  variables 
Q3 .  Qtji  Qg.  0^.  ctc   as  these  have  been  defined  above  in 
Eq.  15  and  in  the  lines  just  following  Eq.  15.   Of  the  158 
reports  listing  the  dependent  sample-data  taken  from  days 
2  through  5  (July  15.  ....  19th,  inclusive)  in  Appendix  B, 
it  was  found  that  eight  cases  had  missing  pyrheliometer 
values  at  the  solar  hour  adjacent  to  the  effective  time  of 
the  NIi'IBUS  scan,  and  thus  led  to  unreasonable  transmissivit- 
ies.   These  pyrheliometer-associated  inconsistencies  in  the 
data  field  were  denoted  by  the  code  symbol  P. 

On  the  other  hand,  for  days  two  through  five,  there  also 
occurred  six  cases  of  negative  QA  after  Eq .  16  was  applied, 
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which  were  also  physically  impossible.   All  cases  for  which 
QA  were  negative  were  coded  by  the  symbol  A  (negative  ab- 
sorptivity).  Hence  the  remaining  sample  of  17^  data  sets 
from  days  two  through  five  were  used  in  the  further  screening- 
analyses  to  be  described  in  Section  III  (C).   For  an  illus- 
tration of  the  initial  screening  involving  data-deletions 
of  code  P  and  A,  reference  should  be  made  to  Table  I  in 
the  column  listing  the  pyrheliometer  station  number.   Sec- 
ondary data-screens  which  lead  to  further  data-exclusions 
based  upon  statistical  tests,  are  utilized  for  additional 
deletions  of  data  cases  based  upon  Improbable  extreme  values 
of 

qA,  qG  and  qR 

These  tests  are  essentially  based  upon  the  requirements: 

qA  «  1  -  qG  (1-Cfc)  -  qR  <^mX)  (l8) 

>qA(MIN) 

Here  the  values  of  q. (MIN)  and  qA(MAX)  are  indicated  by 
empirical  evidence  to  be  discussed  in  Section  III  (C). 

C.   APPLICATION  OF  A  DATA-LIMIT  SCREEN 
1.   General  Considerations 

Hanson  et  al.  (1967)  developed  a  best-fit  regression 
for  q.  using  data  based  upon  TIROS  reflectances,  and  an  em- 
pirically-determined measure  Ify  of  the  water  vapor  absorber 
mass.   Hanson  found  the  most  effective  definition  of  h^f  which 

led  to  the  highest  linear  correlation  of  q.  with  &.,    to  be: 

A       w 

3^ 


Mw  =  (U  sec  Z)*  In  (U  sec  Z)  (19) 

The  specific  form  of  Hanson's  regression  was: 

qA  =  0.096  +  0.045  Mw  (20) 

which  in  turn  was  optimal  in  yielding  the  highest  value  of 
the  correlation  coefficient  R  between  q,  and  My 

R  (q..  MT.)  =  Maximum 

Here  My  of  (19)  was  one  of  many  possible  functions  of  U  sec 
Z  having  the  empirical  form  of  an  absorptivity  (see  Fig.  8). 

Application  of  the  BIMSD  02R  statistical  program  (on 
file  in  the  Program  Library  of  the  W.  R.  Church  Computer 
Center  of  the  Naval  Postgraduate  School)  to  the  dependent 
sample  of  17^  data-cases  from  days  two  through  five,  was 
made  in  order  to  generate  the  least-squares  best-fit  between 

(a)  qA  and  Mw 

(b)  qG  and  Mw 

(c)  qR  and  Mw 

The  most  useful  Immediate  results  obtained  here  were 
listings  of  sample  means  and  standard  deviations  (S.D.)  of 
each  variable  under  consideration,  (as  well  as  listings 
of  the  17^  values)  of  q.,  qc  and  qR  (see  Table  II). 


35 


CO  0.0  10 


Optical  Pathlength,  In  U  Sec  Z 


Figure  8.   Absorptivity  of  solar  radiation  in  the  atmos- 
phere as  a  function  of  optical  pathlength,  In  U  Sec  Z. 
Curve  (1)  is  the  absorptivity  in  a  clear  -atmosphere  by 
H2°»  co2»  and  °2«  fron*  Yamamoto  (1962).   Curve  (2)  is 
the  absorptivity  in  the  clear  atmosphere  by  H2O  only, 
from  Houghton  (195/+)  .   Curve  (3>  is  the  absorptivity  in 
the  atmosphere  over  the  United  States  with  mixed  cloud- 
iness, based  on  the  Hanson's  1967  study.   Vertical  bars 
indicate  the  most  probable  error  at  the  abscissa  value. 
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TABLE  II.   Statistics  Pertaining  to  qA,  qQ,  q 
Mean      S.D.  Range 

qA  0.2290  0.09^3  0.0026  to  0.6192 
qG  O.6308  0.1964  0.07^7  to  0.959^ 
qR    0.2388    0.1127    0.0916  to  0.6406  =  qR(i-iAX) 

2.  Consideration  of  Data-Exclusions  on  q^ 

The  listed  range  of  the  planetary  albedo  q^  in 
Appendix  B  is  within  the  range  of  albedo  values  quoted  in 
the  literature  (List,  1958)  for  various  surfaces:   cloud- 
tops,  soil-types,  etc.   Hence  no  limiting  range  on  the  values 
of  q~  found  here  were  imposed  on  any  of  the  data  cases  of 
this  study,  i.e.,  no  data-case  exclusions  were  made  on  the 

basis  of  qD.   This  does  not  mean  that  q_,  does  not  have  a 
ii  n 

measurement  error,  but  that  its  effect  on  upsetting  the 
balance  of  Eq .  17  is  not  as  great  as  the  more  unrealistic 
ranges  of  values  associated  with  q,  and  q~  . 

3.  Data-Exclusions  Based  on  q. 

Because  qA  was  determined  as  a  residual  involving 
the  two  measured  (but  error-subject)  quantities  qR  and  qG , 
it  was  felt  that  some  reasonable  bounds  would  have  to  be 
applied  to  the  range  of  q,.   Thus  qA(HIN)  and  qA(;iAX)  values 
were  established  in  accordance  with  1%   extreme-probability 
limits  on  both  winces  of  a  Gaussian  distribution  centered  on 
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qA(K£AN).   From  the  statistics  of  Table  II,  It  follows  that 
the  definition  of: 

qA(MIN)  =  0.2290  -  2.3267  (0.0943)  =  0.0096  (21) 

qA(MAX)  =  0.2290  +  2.3267  (0.0943)  ■  0.4484  (22) 

The  value  of  qA(?iIN)  was  relaxed  to: 

qA(MIN)  =0.07  (23) 

[cf.t  Fig.  8,  after  Hanson  et  al.  (1967)]  .   With  the  limits 
established  by  £qs.  22  and  23,  there  were  six  case-deletions 
due  to  small  absorptivity-calculations 

qA<  0.07 

and  4  case  deletions  due  to  unacceptably  large  q. 

qA>  0.4484 

These  ten  deletions  are  coded  A  MIN,  A  MAX  in  Appendix 
Table  B.   The  codes  MIN  and  MAX  indicate  that  statistical- 
empirical  considerations  were  involved  in  these  qA  secondary- 
deletions. 

4.   Data-Exclusions  based  on  q^ 

Here  two  statistical  relationships  analogous  to  Eqs. 
21,  22  for  considering  probable  limits  on  qG  were  set  up: 

qG(MIN)  =  0.6308  -  2.576  (0.1964)  =  0.1249        (24) 

qr(MAX)  =  0.6308  +  2.576  (0.1964)  =  1.1367        (25) 
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The  factors  2.576  In  Eqs.  2*4-,  25  correspond  to  limiting  data- 
cases  lying  In  a  0.5%   rather  than  In  a  lfi   Gaussian  probabil- 
ity distribution  wing.   Here  we  are  being  slightly  sore 
cautious  about  rejecting  extreme  values  of  qG  than  In  the 
corresponding  consideration  of  q„.   Finally  Eq .  25  has  been 
altered  from  an  extreme  transmlsslvlty  of  114^  to  93# 
(I.e.,  0.93) •   This  follows  In  accordance  with  the  minimum 
qA  =  0.07  of  Eqs.  23  and  18  so  that: 

qG(MAX)  =  0.93  (26) 

the  upper  limit  exclusion  £q.  26  accounted  for  no  data- 
case  rejections,  while  the  lower  limit  of  Eq .  2k   accounted 
for  four  data-case  rejections.   These  cases  are  listed  In 
Appendix  Table  B  by  the  symbolism  T  MIN.   The  revised  data- 
sample  size  for  regression  purposes  now  comprises  160  cases. 

Actually  the  purpose  of  reducing  the  number  of  extreme 
data  cases  was  to  eliminate  as  much  noise  as  possible  from 
the  regression  analysis  so  that  the  physical-statistical 
relationships  would  stand  out  more  clearly.   It  will  be 
seen  in  Section  IV  that  the  data-screening  process  for 
elimination  of  extreme-cases  did  not  improve  the  regression 
analysis  markedly,  after  the  initial  reduction  to  the  sample 
of  174  cases  (based  upon  obviously  missing  and  erroneous 
data)  was  achieved. 
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IV.   THE  REGRESSION  RESULTS 

A.   THE  REGRESSION  ANALYSIS  FOR  qA 

Initially  a  least-squares  fit  analogous  to  that  of 
Hanson  (Eqs.  19.  20)  was  sought,  using  both  the  full  sample- 
set  of  17^  cases,  and  the  smaller  sample-set  of  160  cases 
after  extreme-value  limits  had  been  applied  to  screen 
statistically  for  possible  unreliable  data.   The  regression 
method  used  was  that  of  3II1ED  02R  (Dixon,  1970).   Seven 
variables  were  listed  on  an  IBM  card  for  each  data-case. 
Following  are  the  variables  listed  for  each  case  in  Appendix 
Table  B,  namely: 

Qs.  QR.  Qq.  Qa,  U,  cos  Z,aG 

By  means  of  transgeneration  (an  option  of  the  BIKED  series 
of  programs),  the  fractions  q,  =  QA/Qs»  ^g •  °*R  were  formed, 
as  well  as  the  water  vapor  mass  M.«  (Eq.  19)  used  by  Hanson 
et  al.  (1967). 

Two  types  of  regression  tests  were  run  for  both  the  raw- 
data  and  screened-data  sample  sets.   These  tests  involved: 

(1)  qA  on  Mw  alone 

(2)  qA  on  both  M^  and  qR 

The  reason  for  performing  the  regression-test  (2)  is  that 
it  was  observed  that  an  above  (below)  normal  qA  tended  to 
occur  with  an  above  (below)  normal  qD.   In  fact  the  correla- 

h 

tion  coefficient  between  qA  and  q  was  R (q . ,  q  )  =  0.335* 
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More  specifically,  an  Increase  In  qu  Is,  In  general,  as- 

n 

soclated  with  an  Increase  of  cloud  cover,  so  that  the  In- 
clusion of  qR  as  a  predictor  for  q,  should  Indicate  at 
least  partially  an  effect  of  the  cloud  cover  on  q  .   riore- 
over  with  the  NOAA  satellite  scanning  radiometer  readout 
available  at  Fleet  Weather  Centrals,  values  of  qR  ought  to 
be  as  readily  accessible  as  the  radiosonde  data  used  for 
the  M.  computation. 

The  two  regression  tests  on  the  raw  and  statistically- 
screened  data  sets  led  to  the  following  regression  equations 
and  associated  correlation  coefficients  R. 

Test  1,  sample  size  n  =  17^  (raw  data  test): 

qA  =  0.1983  +  0.0188  My  (2?) 

R  =  0.1807 

Test  1,  sample  size  n  =  160  (screened  data): 

qA  =  0.1978  +  0.0189  Mw  (28) 

R  =  0.2206 

Test  2,  sample  size  n  =  17^  (raw  data): 

qA  =  0.16*19  +  0.1880  qR  +  0.0118  KM  (29) 

R  =  0.2802 

Test  2,  sample  size  n  =  160  (screened  data): 

qA  =  0.1605  +  0.2162  qR  +  0.0111  rtrt  (30) 

R  =  0.3567 


*U 


2 

Since  the  correlation  coefficient  squared  (R  )  is  the 

measure  of  the  explained  fraction  of  the  variance  Intro- 
duced by  the  regression  technique,  it  is  clear  that  q*  Is 
not  well  specified  by  either  St.   and/or  by  q^  in  spite  of 
evidence  to  the  contrary  by  Hanson  et  al.  (1967).   The 
most  significant  correlation  for  q-  of  the  set  of  Lqs. 
27,  28,  29,  30,  is  that  involving  both  q^  and  H^  as  pre- 
dictands  with  statistically  screened  data.   The  mean  sta- 
tistics which  relate  to  the  four  tests  Just  cited  are 
summarized  in  Table  III.   The  mean  standard  error  of  es- 
timate for  the  best  case  (line  k)    in  Table  III  was  S.S.  = 
0.0715. 

The  relatively  large  standard  error  which  remains  after 
the  regression  tests  (Table  III),  seems  to  be  due  to  the 
computational  formula  for  q.  of  (17) t 

*A  =  1   -  qR  -  (1-aG^G  (17) 

This  equation  has  been  used  with  partly  cloudy  to  cloudy 
skies,  when  qG  has  already  been  reduced  by  the  cloudiness 
associated  with  the  large  reflectance  q^.   However,  in 
these  cloudy  cases,  the  multiplicative  reduction  of  the  qG 
term  of  (17)  by  the  factor  (1-a,-)  gives  a  positive  bias 

u 

to  the  computation  of  q»  which  shows  up  in  the  predominance 

A 

of  large  positive  residuals  remaining  after  application  of 
the  regressions,  Eqs.  27,  28,  29,  30. 
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Table  III 

Regression  Statistics  for 
Tests  of  (1)  qA  on  Mw;  (2)  qA  on  My,  qR 

Standard 

Sample  Mean  Deviation  Corr.  %   ExpL  Std.  Err  . 

Test        Size    Qa    °  IgAJ    Coeff.  feiiance  of  Est. 

(1)  qA  on  M^      174  0.2290  O.09U3  0.1807  3-27  0.0931 

(1)  qA  on  Mw      160   0.2275   0.0760   0.2206  4.87   0.0744 

(2)  qA  on  %,qR   174   0.2290   0.0943   0.2802  7.85   0.0911 
(2)  qA  on  Mw,qR   160  0.2275  O.076O   0.3567  12.72   0.0715 
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B.   A  REGRESSION  ANALYSIS  FOR  qG 

The  relatively  weak  statistical  specification  of  q   was 
due  In  part  to  Improper  accounting  for  earth-surface  absorp- 
tion term  In  Eq.  17  when  clouds  were  present.   However  qG 
Is  a  measured  parameter,  which  Is  related  to  q.  through 

Clear  case:    qG  =  i'^R"^  (3D 

^G 

Overcast  case:  qG  =  l-qR-q,  (32) 

Hence,  since  the  computational  error  In  determining  qA  tends 
to  obscure  the  statistical  regression  of  qA  on  My  and  q., . 
Eqs.  31i  32  suggest  that  such  a  statistical  relationship 
might  be  contained  In  a  regression  of  the  transmlsslvlty 
qr  on  M.r  and  q„  In  view  of  the  functional  relationships  (3D, 
(32)  between  q~  and  q,.   There  seems  to  be  a  greater  likeli- 

u        A 

hood  of  a  significant  regression  for  qr  than  for  q   in  view 
of  the  fact  that  q   is  directly  measured.   The  assumption  of 
a  significant  regression  of  the  type  just  described  is  the 
hypothesis  of  subsection  IV  (B). 

In  the  statistical  tests  to  be  conducted  here,  qG  is 
established  as  the  predictand  in  a  set  of  regression  tests 
for 

(1)  qG  on  Mu  alone, 

(2)  qG  on  My  and  q^. 

These  tests  then  led  to  the  following  regression  results  in 
Eqs.  33,  34,  35,  36  which  are  the  analogs  of  Eqs.  27,  28,  29, 
30,  respectively.   The  associated  regression  statistics  are 
listed  In  Table  IV. 
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Table  IV 

Regression  Statistics  for 
Tests  of  (1)  qG  on  Mw;  (2)  qG  on  My  and  qR 

Std. 
Sample  Mean    9evr    Corr.  %   Sxpl.  Std.  Err. 
Test       Size    qG     ^Up)    Coeff .  Variance  of  Est. 

(1)  qG  on  Mw  174  0.6308  0.1964  0.3124  9.77  0.1871 

(1)  qG  on  Mw  160  0.6429  0.1752  0.3369  11-35  O.1654 

(2)  qG  on  I%.qR  174  O.6308  0.1964  0.6236  67.83  0.1120 
(2)  qG  on  Mw.qR  160  0.6429  0.1752  0.8679  75-32  O.O876 
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Test    1,    sample   size   n  =   174    (raw  data): 

qG   =   0.74l4l   -   0.0677/irt  (33) 

R  =   0.3124 

Test   1,    sample   size   n  =   160    (screened  data): 

qG   =   0.747^3   -   0.066514^  (34) 

R  =  0.3369 

Test  2,  sample  size  n  =  174  (raw  data): 

qG  =  0.9887  -  1.3927qH  -  0.0155^  (35) 

R  =  0.8236 

Test  2,  sample  size  n  =  160  (screened  data): 

qG  =  0.99252  -  1.420?8qR  -  0.01509.^  (36) 

R  =  0.8679 

Note  the  relatively  high  correlation  coefficients  in 
the  Test  (2)  case  in  the  specification  of  qG  in  terms  of 
both  I4W  and  qR .   The  results  seem  to  be  slightly  better  in 
the  case  of  filtered  data  (sample  size  of  160)  ,  where  the 
standard  error  of  estimate  of  qG  is  O.O876  compared  to 
0.1120  in  the  case  of  the  raw  data. 

In  Table  IV,  the  lower  values  of  the  correlation  of  qG 
on  LK  alone  suggests  that  the  qG  estimations  developed  from 
(33).  (34)  will  not  be  as  useful  as  those  derived  from  (35). 
(36).   This  will  be  further  explained  in  Subsection  (C ) . 

C.   TESTS  ON  INDEPENDENT  DATA -SAMPLES 

In  Sections  IV  (A)  and  (b) ,  four  sets  of  regression 

equations  for  both  q;,  and  q  were  statistically  deduced 

•^       G 
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based  upon  the  dependent -data  samples  of  17^  and  160  cases, 
respectively  (termed  the  raw-data  and  statistically-screened 
samples,  respectively). 

The  regression  equations  so  derived  were  of  the  form: 

*qA  =  AQ  +  A1qH  +  A^  (37) 

-qG  =  BQ  +  BxqR  +  Bgl^  (38) 

and  represent  approximations  to  the  observed  or  measured 
q»,  qG.   In  Eqs.  27  and  26,  the  special  case 

A-l  =  0 

was  enforced,  while  in  Eqs.  33 t  3^.  we  required 

B1   =  0 

Eqs.  37t  38  with  A,  =  0,  B  =  0,  respectively  were  termed 
the  Test  1  regressions,  while  A,  ^  0,  B,  ^  0,  have  been 
termed  the  Test  2  regressions. 

The  generating  formulas  of  Eqs.  37  and  38  may  be  written 
in  the  matrix-product  notation: 


vV&) 


^qA  =  (A0,  Av  A2)  [qqj  (39) 


B 


q   =  (B0.  B  ,  B2)  [qR\  (40) 


Here  AQ,  A^  a2  are  listed  in  Eqs.  37,  38,  39.  ^0  and  B  , 
Bi,  B2  are  listed  In  (33 36)". 


k? 


Tabl( 

5  V 

Independent  Sample 

Regression 

Tests 

(a)   Raw-Data  Sample 

1.  n=44 

Predic- 

tand 

Heans 

Standard 
Dev.a(q) 

Standard 
Error 

Correlation 
Coefficient 

qA=0.2179 

0.0882 

qA 

lfA=0.228l 

0.0263 

0.0860 

R=0.2652, 

by 

Eq  . 

27 

^=0.2268 

0.0290 

0.0871 

R=0.2155, 

by 

Eq  . 

29 

qG 

qG=0.6575 

0.1718 



fG=0.6372 

0.0805 

0.1688 

R=0.2387. 

by 

Eq. 

33 

^=0.6465 

0.1544 

0.1073 

R=0.7868, 

by 

Eq. 

35 

(b) 

Screened  Data 
qA=0.2124 
fA=0.2253 

.,  n=4l 
0.0688 
0.0198 

0.0688 

qA 

R=0.0551, 

by 

Eq. 

28 

qA=0.2240 

0.0246 

0.0688 

R=0.0975. 

by 

Eq. 

30 

qG 

qG=0.6745 

0.1467 

^G=0.6508 

O.O696 

0.1467 

R=0.0665, 

by 

Eq . 

34 

^=0.6588 

0.1383 

0.0882 

R=0.8048, 

by 

Eq. 

36 
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For  the  Independent  test,  the  simultaneous  samples 
(qA,  q  )  and  (qG ,  *q\)  from  day  1  were  used.   For  the  test 
on  the  raw -data  Independent  sample  of  day  one,  the  four 
specific  coefficient  sets  of  A_,  A.,  A2  and  5q,  B.  ,  B£  of 
Eqs.  27,  29,  33,  35  were  used.   For  the  test  on  the  screened- 
data  Independent  sample  of  day  one,  the  four  coefficient 
sets  of  Eqs.  28,  30,  3^,  36,  respectively  were  used.   Ihe 
regression  results  for  the  eight  combinations  of  the  one- 
predictor  and  two-predictor  formulas  for  q,  and  q,.  are  shown 
in  Table  V  for  both  the  raw  and  screened  independent  sample 
tests. 

Table  V  shows  that  the  use  of  screened  data  in  place  of 
raw  data  caused  a  decrease  In  the  standard  deviation  of  the 
predlctand  by  eliminating  some  of  the  extreme  oalues.   r.ore- 
over  the  correlation  coefficients: 

R(qA.  qA).  R(qG.  %) 

decreased  from  the  raw  to  the  screened  cases  when  the  listed 
values  of  the  correlation  coefficients  of  the  initial  spec- 
ification Eqs.  27,  28,  29,  30,  33,  3*4-  were  small  (in  the 
range  .2  to  .35)*   However  the  two  key  results  which  stand 
out  in  Table  V  are  those  corresponding  to  the  use  of  specify- 
ing Eqs.  35,  36  (which  gave  R  =  0.7868  and  0.8048,  respec- 
tively), denoted  by  an  asterisk  superscript  in  Table  V. 

These  two  key  results  indicate  that  the  transmissivlty- 
estimator  "q^  expressed  by  either  Eq .  35  or  36  gives  a  sig- 
nificantly good  specification  for  the  pyrheliometric  obser- 
vation of  qG. 
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The  Interpretation  of  the  correlation  coefficients  R 
In  the  two  cases  noted  In  Table  V  may  be  seen  from  the 
standard,  statistical  formula: 

R2  _  Mean  square  explained  by  Regression  (Ul) 

Mean  square  of  predlctand 

aq      -  (SE)2 

Thus  by  use  of  Eq.  35 .  the  standard  error  is  0.1073t 
whereas  the  original  standard  deviation  is  0.1718.   Like- 
wise by  use  of  £q .  J6   the  standard  error  is  0.0882  compared 
to  the  original  standard  deviation  of  0.146?,  as  summarized 
by  Table  V. 

Actually,  there  has  been  no  appreciable  significant 
change  in  the  fractional  part  of  explained  variance  by  (35) 
versus  (36).      Moreover  the  process  of  providing  screening 
limits  is  very  laborious  if  it  provides  no  significant  im- 
provement in  the  fractional  explained  variance  (right  side 
of  Eq .  41 ) . 

D.   ERROR  CONSIDERATIONS 

As  has  already  been  noted,  there  Is  a  substantial  stand- 
ard deviation  in  q,„  (0.1964  in  the  case  of  raw  data,  0.1752 
in  the  case  of  screened).   A  substantial  part  of  this  error 
is  due  to  initial  calibration  of  the  Eppley  pyrheliometers , 
ranging  from  3*5  to  10. 5^ t  Hanson  et  al.  (1972).   While 
such  errors  may  be  partially  accounted  for  and  correction 
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factors  applied,  the  nature  of  the  correction  process  was 
not  known  until  the  present  study  was  completed.   It  is 
likely  that  a  sizeable  fraction  of  the  standard  error  of 
estimate  of  qG  listed  in  the  last  line  of  Table  IV  is  due 
to  the  pyrheliometer  calibration  error. 

Another  minor  source  of  instrumental  error  is  that  due 
to  the  necessity  of  interpolation  of  qG  in  time  between 
tabulated  hourly  solar  transmission  increments.   The  non- 
availability of  detailed  digitized  records  of  the  pyrhelio- 
meter and  satellite  reflectance  data  both  on  a  sensitive 
time  scale  was  another  error-source  in  this  study  since 
short-term  cloud  variations  of  a  period  less  than  a  half- 
hour  would  have  been  related  to  the  short  term  variabilities 
of  qG  and  of  qH.   Note  that  in  this  connection,  our  time- 
smoothed  records  of  qG  and  q^  were  highly  correlated,  there 
being  a  correlation  coefficient  R  (qG,  qfi)  =  -0.821. 

As  mentioned  earlier  in  connection  with  the  effect  of 
cloud-cover  on  the  computation  of  qG  (Eqs.  31 ■  32),  two 
different  formulas  apply  in  the  clear  and  overcast  cases 
for  the  solar-insolation  budget  as  accounted  for  at  the 
pyrheliometer.   Failure  to  explicitly  account  for  the  amount 
of  cloud-cover  has  biased  the  computation  of  q^  in  the  sense 
of  increasing  the  value  of  the  latter  by  up  to  15/-  (that 
is,  by  the  factor  CL/(l-aG)  In  Eq .  31,  where  CL  is  the  cloud- 
cover).   A  discussion  for  the  systematic  inclusion  of  the 
cloud  cover  parameter  for  improved  specification  of  qA  and 
qG  is  given  in  Section  V. 
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V.   RECO/IHENDATIONS  AND  CONCLUSIONS 

As  deduced  from  this  study,  the  insolational  budgetary 
model  of  Eqs.  17  and  31  gave  q.  values  biased  towards  high 
values.   This  was  due  to  the  failure  to  account  for  the 
effects  of  cloud  cover  properly,  as  shown  for  the  overcast 
case  In  Eq.  32.   The  combination  of  Eqs.  31 .  32  show  that 
for  some  intermediate  cloud-cover  state  CL,  the  properly 
weighted  form  of  these  two  equations  leads  to 

qA  =  1   -   qR  "  qG  t1  "aG  t1-01^]  =  °  (i+2> 

In  Eq.  42,  the  term  Involving  a^  is  the  fractional  insola- 
tion which  the  cloud-openings  (1-CL)  permit  to  escape  to 
space. 

In  this  study,  CL  was  not  known  and  was  assumed  zero, 
as  was  also  done  by  Hanson  et  al.  (1967)  and  Fritz  et  al. 
(1964).   However,  recent  studies  by  Downey  et  al.  (1972) 
showed  that  qR  is  highly  correlated  with  CL,  as  indicated 
by  the  following  summary  of  his  results  (based  on  NIMBUS  III 
observations) . 

Correlation  coefficients  R(qR,  CL)  by  periods  1969-70. 

May  16-31    Aug.  1-15    Oct.  3-17    21  Jan-3  Feb. 

_____    '  1970 

o7o"5       0.78       0.78        0.85 

With  results  of  this  type,  combined  with  satellite 
measurements  of  window-channel  temperatures  as  well,  seasonal 
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regressions  of  some  physically  realistic  form  should  be 
fitted  to  the  observed  cloud-coverage.   The  objective  is 
to  obtain  a  highly  specifying  representation  for  CL  of  form 

CL  =  F(qR.  TBfi)  (43) 

where  F  is  the  desired  statistically-generated  function  of 
the  indicated  variables. 

A  seasonal  variation  may  be  built  into  the  function  F 
by  making  F  dependent  upon  sec  Z,  where  Z  =  zenith  angle. 

Finally  the  observed  values  of  q,~  which  are  used  in  con- 
nection with  the  solution  of  (k2)  for  qA  should  be  carefully 
corrected  for  initial  calibration  and  pyrheliometer  degrada- 
tion errors  after  the  procedure  recommended  by  Hanson  et  al. 
(1972). 
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APPENDIX  TABLE  A 
NORTH  AMERICAN  PYRHELIOMETER  STATION  INFORMATION 


Stat. 

No. 

Code 

1 

AA 

2 

AB 

3 

AC 

4 

AD 

5 

AE 

6 

AF 

7 

AG 

8 

AH 

9 

AI 

10 

AJ 

RS 

11 

AK 

12 

AL 

13 

AM 

14 

AN 

RS 

15 

AO 

RS 

16 

AP 

RS 

17 

AQ 

18 

AR 

19 

AS 

20 

AT 

RS 

21 

AU 

RS 

22 

AV 

23 

AW 

24 

AX 

RS 

25 

AY 

26 

AZ 

Name  Index 

Fort  Worth,  Texas  72259 

Apalachlcola.  Fla.  220 

Miami,  Florida  202 

Brownsville,  Texas  250 

Charleston,  S.C.  208 

Nashville,  Tenn.  327 

Columbia,  Mo.  445 

Dodge  City,  Kansas  451 

Caribou,  Maine  712 

Madison,  Wise.  64l 

Green  Bay,  Wis.  645 

El  Paso,  Texas  2?0 

Albuquerque,  N.M.  365 

Ely,  Nevada  486 

Los  Angeles,  Ca.  295 

San  Diego,  Ca.  290 

Phoenix,  Ariz.  278 

Tucson,  Ariz.  274 

Santa  Maria,  Ca.  394 

Vandenburg,  Ca.  393 

Bismark,  N.D.  764 

Great  Falls,  Mont.  775 

Medford,  Oregon  597 
Tacoma,  Wash. 

Tatoosh,  Wash.  798 

Fresno,  Ca.  389 

Oakland,  Ca.  493 

Cape  Hatteras,  N.C.  304 
Sterling,  Va. 

Boston,  Mass.  509 

Nantucket,  Mass.  506 

New  York,  N.Y.  74486 

Omaha,  Nebr.  72553 


Lat .   Long.  Albedo 


32-50 
29-44 
25-48 
25-54 
32-54 
36-07 
36-58 
37-46 
46-52 
43-08 

31-48 
35-03 
39-17 
33-56 


46-46 

47-29 
42-22 

47-57 


35-16 
39-00 
42-22 

40-39 
41-22 


97-03 
84-59 
80-16 
97-26 
80-02 
86-41 
92-22 
99-58 
68-01 
89-20 

106-24 

106-37 
114-51 

118-23 


33-26     112-01 


100-45 
111-22 
122-52 
124-33 


75-33 
77-26 
71-02 

73-47 
96-03 


16 
12 

15 
16 

13 
14 

16 
18 
14 
17 

19 
22 

19 
15 

22 


34-54   120-27    15 


15 
17 
16 
16 


36-46  119-43   13 


14 
15 
15 

15 

15 
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Stat 

• 

No. 

Code 

27 

CA 

RS 

28 

CB 

29 

CC 

RS 

30 

CD 

31 

CE 

RS 

32 

CF 

33 

CG 

RS 

34 

CH 

RS 

35 

CI 

RS 

36 

CJ 

RS 

37 

CK 

RS 

38 

CL 

RS 

39 

CM 

RS 

4o 

CN 

RS 

41 

CO 

RS 

42 

CP 

RS 

43 

CQ 

RS 

44 

CR 

RS 

45 

CS 

RS 

46 

CT 

RS 

47 

CU 

Name 


Index   Lat .   Lon*r.  Albedo 


RS 


Beaverlodge  CDA,  Alta. 

55-12 

119-30 

14 

Prince  George  b.C. 

72896 

Churchill,  Man. 

913 

58-45 

90-05 

15 

Departure  Bay,  B.C. 

49-12 

124-00 

15 

Iatoosh,  Wash. 

798 

Edmonton,  Alta, 

74119 

53-19 

113-35 

15 

Fredericton  CDA,  N.B. 

45-52 

66-32 

15 

Caribou,  ME 

72712 

Goose  A,  (Nfld) 

816 

53-19 

60-25 

16 

Guelph  OAC  (Ont. ) 

43-36 

60-18 

14 

Buffalo,  N.Y. 

526 

Halifax,  N.S. 

44-53 

63-31 

16 

Sable,  ISC 

600 

Kapuskaslng  A.  Ont. 

49-25 

82-28 

15 

Moosonee  CDA,  Ont. 

836 

Kentville  CDA  N.S. 

45-06 

64-30 

15 

Portland,  ME 

606 

Montreal,  Qua. 

45-30 

73-42 

16 

Manlwakl,  CDA,  Que 

722 

Norinand  in,  C DA  ,  Que 

4G-54 

72-36 

15 

Sept.  Isles,  CDA,  Que 

811 

Ottawa  NRC ,  Ont . 

45-24 

75-42 

16 

Manlwakl,  CDA,  Que 

722 

St.  John's  West, CDA, (Nfld] 

1 

47-37 

52-45 

16 

Stephenvllle,  (Nfld) 

815 

Suffield  A,  Alta. 

50-16 

111-11 

15 

Great  Falls,  Mont. 

775 

Summer land  CDA,  B.C. 

49-36 

119-42 

16 

Spokane,  Wash. 

785 

Swift  Current  CDA,  Sask. 

50-17 

107-41 

15 

Glasgow,  Mont. 

768 

Toronto,  Ont. 

43-43 

79-14 

14 

Buffalo,  N.Y. 

528 

Toronto  Scarborough,  Ont. 

43-43 

79-14 

14 

Buffalo,  N.Y. 

528 

Vancover  UBC ,  B.C. 

49-11 

123-10 

16 

Tatoosh,  Wash. 

798 

Winnipeg  Int  A.  Man. 

49-54 

97-14 

13 

International  Falls, Minn. 

747 
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APPENDIX  B 

Tables  of  observational  data  Interpolated  to  the  effec- 
tive scan-time  at  each  of  the  k7   stations  used  in  this  study, 
and  for  days  1,  ....  5.  (July  15 19.  respectively). 

The  symbols  Qs,  QH,  G~  ,  Q.,  a^   are  as  defined  in  £q.  15 
and  in  the  paragraph  which  Immediately  follows  this  equa- 
tion.  The  symbol  U  denotes  the  precipitable  water  vapor 
at  the  station  and  CZ  is  the  cosine  of  the  zenith  angle  at 
the  effective  time  of  the  observation. 

DELETE  CODS  symbols  P  and  A  Indicate  inconsistent  trans- 
mission and  absorption  values,  respectively,  and  these  data- 
cases  were  not  used.   All  other  codes,  such  as  A  i'^AX,  a  KIN. 
T  MIN  express  the  occurrence  cf  statistically  Improbable 
values,  which  cases  were  used  in  a  secondary  data-cases 
(Sec.  III.  (C)). 
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